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ABSTRACT 


The  Fisher  detector  and  conventional  power  detector  have 
been  implemented  -n  short-period  data  recorded  at  .he  Korean  Seismic  Re- 
search Station.  It  was  found  that  false  alarm  rates  varied  so  widely  from 
day  to  day  that  no  comparison  between  detectors  could  be  made.  However, 
a comparison  between  detectors  with  different  integration  times  but  the  same 
false  alarm  rate  revealed  ihat  integration  time  does  not  affect  the  probability 
of  detection.  A means  to  hold  the  alarm  rate  constant  should  be  developed  to 
permit  comparison  and  evaluation  of  these  detectors. 


Neither  the  Advanced  Research  Projects  Agency  nor  the  Air  Force 
lechnical  Applications  Center  will  be  responsible  for  information  contained 
herein  which  has  been  supplied  by  other  organizations  or  contractors,  and  this 
document  is  subject  to  Later  revision  as  may  be  necessary.  The  view's  and  con- 
clusions presented  are  those  of  the  authors  and  should  not  be  interpreted  as 
necessarily  representing  the  official  policies,  either  expressed  or  implied,  of 
the  Advanced  Research  Projects  Agency,  the  Air  Force*  Technical  Applications 
Center,  or  the  US  Government. 
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1 his  repo rt  presents  results  of  a study  of  two  seismic 
signal  detectors;  the  Fisher  detector  and  the  conventional  power  detector. 
These  detectors  have  been  described  in  an  earlier  report  (Lane,  1973', 
herein  called  Report  1,  where  their  performance  on  long-period  data  was 
studied.  The  present  report  is  concerned  with  their  response  to  short- 
period  data  at  the  Korean  Seismic  Research  Station  (KSRS1. 

A total  of  185  events  were  processed  by  the  detectors  and 
the  results  used  to  estimate  a probability  of  detection  for  magnitudes  greater 
than  3.0.  A number  ot  hour-long  noise  samples  were  also  processed  to  find 
false  alarm  rates.  A simple  qua  1 ity -cont  rol  algorithm  was  devised  to  re- 
move spiked  and  dead  channels. 


SECTION  II 


ALGORITHMS,  DEFINITIONS,  AND  DATA 
A.  ALGORITHMS 

Roth  the  Fisher  and  conventional  power  detector  reported 
on  here  have  been  described  from  a theoretical  and  experimental  point  of 
view  in  Report  I.  Each  searches  the  time-delayed  data  for  signal-like 
characteristics.  In  the  case  of  the  Fisher  detector  this  characteristic  is 
waveform  similarity.  If  the  time-aligned  traces  are  designated  by 
y.,  the  Fisher  detector  output  is 
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where  M is  the  number  of  channels  and  the  bar  denotes  an  equally  weighted 
average  over  some  previous  time  interval,  called  the  integration  time. 

When  all  the  y(  are  the  same,  the  denominator  of  II- 1 is  zero,  because  under 
these  conditions  the  mean  square  of  the  y.  is  equal  to  their  squared  mean. 
When  the  y.  are  perfectly  uncorrelated,  the  detector  output,  as  normalized, 


in  line.  A real  signal  will  al»a\  s be  accompanied  b^  noise,  and  for  the 
rase  of  propagating  undistorted  plane  wave  signals  with  noise  correlated 
neither  with  itself  nor  the  signal,  the  detector  output  can  be  shown  to  be 
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■here  S is  the  signal  (lower  .end  N is  the  noise  power. 


Rea!  signals  .ire  generallv  dissimilar  from  site  to  site  .is 
‘•'II  as  noisy,  and  the  Fi  sher  detector  lUtput  is  reduced  in  this  ease- 
i.xpr sooioi.s  lor  this  reduced  output  are  derived  in  Report  l,  where  it 
is  sh>v.  n that  the  detector  saturates  at  an  output  le-vel  dependent  on 
t ■'  «•  hind  are!  degree  ,»i  distortion. 


1 he  e line  ntlem.ll  power  dc  tec  tor  to  rills  the  ratio  of  the- 
-hort-te  rm  av.  rage  beam  (lower  to  the  long-term  average  beam  powe  r 
Ihe-  short-term  average  is  equally  weighted  aveT  its  time  gate,  while-  th< 
long-te  rm  ave  rage  is  exponent  iall  v weighted  into  the  past.  Ihe  detector 
• I 'll  put  is  t lien 
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h«  re  IT  is  the  I mg-te  rm  .ever  ege  anei  the  otlier  te  rms  are  defined  as  abov 
1-or  undi  storteel  plane-  wave  propagation  the  output  is  again  equal  to  equation 
l!-d  an  el  the  limiting  values  for  zero  signal  an<l  zero  noise  .ere  the-  same.  II 
« oiiv*  n tiona  I dete-clu  r is  li  ‘.s  allee  t e -el  hv  ele- vi.’i  t ion  trom  plane  wave  liehavior 


than  1 he  Fisher  detector,  however,  since  it  merely  require  s that  the  power 
i n < re-use-  on  .ill  channels,  ruthc-r  than  that  the  channels  he-  similar.  It 
saturates  at  a highe-r  level  for  the-  same  degree  and  kinel  of  elistortion. 

Both  tlie-  Fisher  and  conventional  power  detectors  involve  a 
number  of  rhn  rae  le-ristie  linie-s.  These-  include*  the  upper  nnel  lowe  r cut- 
oil  pe-riods  of  the-  pre  filter,  the-  long-term  averaging  time  for  the  conven- 
tiemal  de-tee  tor,  and  llie-  inte  gration  times  for  the-  elete-ctor  outputs.  To 
I >i  e k tentative  values  for  ilu-se-  time-s,  each  time  was  scaled  from  its  value 
in  the-  lone-period  detector  by  the  ratio  of  the  pe-riods  of  pe-ak  motion  in 
the-  .short  and  lone-pe  riod  data.  This  ratio  was  about  .02S,  so  starting 
values  for  the-  intc-eration  times  were  chosen  to  be  0.8  to  6.  -4  seconds, 
and  the  data  were  band-pass  filte-re-el  be-twe-en  0.  1 and  2.0  se-conels.  The 
time  constant  for  the-  long-term  ave-rage  was  reduced  to  6 se-conels. 

A simple  epial  i ty  cheek  was  devised  to  remove  faulty  data. 

In  this  algorithm  the  power  for  each  channel  was  average-el  over  about  2d 
se  conds,  and  ihe-  median  e banned  power  found.  Channels  with  powers 
diffe  ring  from  this  value-  by  more-  than  a factor  of  1 in  e-ithe-r  direction  were 
discarded.  This  procedure-  was  as  successful  at  detec  ting  spike-el  ox-  de-ad 
c hannels  as  inspection  of  the-  time-  traces  in  all  ease-s  whe-i’e-  the  traces 
w e re  examine  el. 

B.  DEFINITIONS 

The  detectors  stuelie-d  he  re-  involve-  only  one-  output,  as  defined 
in  equations  II - 1 nnel  II- 1.  When  this  output  rises  above  some-  pr e-determined 
value,  called  the-  detec  tion  thre  shold,  it  is  claimed  that  a seismic  signal  is 
pre  sent.  A log  of  such  claimed  detections  is  ke-pt,  and  an  analyst  may  later 
confirm  or  re-jeet  each  one-,  although  this  was  not  done  in  the  pre-sent  report. 
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The  detector  output  may  rise  above  the  threshold  for  some 
reason  other  than  the  presence  of  a real  signal.  This  is  called  a false 
alarm,  and  an  assessment  of  the  rate  at  which  these  alarms  occur  is  an 
important  part  of  the  detector  evaluation. 

Conversely,  the  detector  may  not  respond  sufficiently  to  a 
real  signal,  due  to  distortion  or  masking  noise.  The  probability  of 
detection  of  a class  of  events  is  the  fraction  of  those  events  for  which  the 
detector  output  rose  above  some  threshold.  This  is  also  an  important 

detector  characteristic  and,  like  the  false  alarm  rate,  is  dependent  on 
threshold  setting. 


C. 


DATA 


The  Korean  Seismic  Research  Station  (KSRS)  is  located  in 
Korea,  at  37°  27'  North  latitude  and  127°  54'  East  longitude.  It  consists 
Of  a 19  element  hexagonal  array  of  short-period  instruments  arranged 
in  two  concentric  rings  about  a center  element. 

Inter-element  spacing  is  about  2.  5 kilometers,  roughly  half 
that  at  the  NORSAR  short-period  array.  KSRS  also  has  an  array  of  seven 
long-period  instruments;  however,  these  instruments  were  not  used  in 
this  study.  The  short-period  data,  for  the  vertical  component  only.are  re- 
corded on  magnetic  tape  at  the  rate  of  20  points  per  second. 

Data  are  available  on  library  tapes  of  8 hour  duration,  and 
one  such  tape  was  present  for  every  alternate  day  starting  29  April  1973, 
and  ending  26  July  1,73.  This  represents  360  hours  of  data  over  a three’- 
month  period.  During  the  time  covered  by  the  tapes  185  seismic  events 
of  magnitude  equal  to  or  greater  than  3.0  and  of  epicentral  distance  less 
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than  80  degrees  but  greater  than  20  degrees  from  KSRS  were  reported  in 
the  NORSAR  and  LASA  bulletins.  These  events  form  the  data  base  for  this 
report.  A histogram  of  their  distribution  with  magnitude  is  shown  in  Figure 
II-l,  and  their  distribution  by  region  is  shown  in  Table  II- 1. 

Events  of  epicentral  distance  less  than  20  degrees  were 
excluded  from  this  study.  At  this  distance,  waveform  similarity  from  site 
to  site  is  poor,  and  time  delays  are  calculated  less  accurately  from  plane 
wave  models  than  for  more  distant  events.  Thus  the  Fisher  detector,  and  to 
a somewhat  lesser  extent  the  conventional  power  detector,  would  not  be 
expected  to  respond  to  these  events. 

These  events  are  not  within  the  area  of  interest  since  they 
are  largely  in  the  Japanese  Islands.  If  they  were  included  in  the  data  base, 
they  would  significantly  bias  the  detection  probability  downward  since  there 
are  about  two  thirds  as  many  events  reported  within  20  degrees  as  beyond. 
Therefore,  these  events  have  been  treated  as  noise.  They  have  not  been 
included  in  the  ds.ta  base  as  events,  nor  excluded  from  the  "signal-free" 
noise  samples.  Their  presence  in  the  noise  samples  may  increase  the 
false  alarm  rate,  but  will  do  so  just  as  in  an  operational  detector. 

Examination  of  the  NORSAR,  and  LASA  bulletins  revealed  that 
no  events  were  reported  during  the  following  time  intervals: 

Day  Hours 

131  0 3:00:00  to  04:00:00 

141  1 3:  30:00  to  14:  30:00 

161  13:  30:00  to  14:  30:00 

183  09:00:00  to  10:00:00 

These  samples  were  used  to  calculate  noise  statistics. 


t 

i 


22  H 


10 

t— 

Z 

41 

> 

UJ 

b 

Qc 

& 

xr 

3 


30  - 
10  - 
16  - 
14  - 
1Z  - 
10  - 
0 - 
6 - 
4 - 


MAGNITUDE  Cue) 


FIGURE  II- 1 

HISTOGRAM  OF  NUMBER  OF  EVENTS  BY  MAGNITUDE 
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The  quality  of  the  data  was  low.  Dead  channels  and  spikes 
were  common,  but  these  data  were  removed  to  a large  extent  by  the  quality 
check.  From  as  few  as  7 to  as  many  as  all  19  channels  were  available  at 
various  times. 

Delay  times  found  by  visually  aligning  traces  agreed  well 
with  those  predicted  by  a plane  wave  model.  However,  amplitude  distortion, 
probably  due  to  the  recording  system,  was  sometimes  noted  and  has  been 
observed  by  others  as  well  (Der,  1974). 
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SECTION  III 


FALSE  ALARM  RATE 


To  find  the  response  of  the  detectors  to  pure  noise,  the  hour- 
long  noise  samples  mentioned  in  Section  II  were  processed  with  each 
detector.  The  detectors  used  integration  times  of  0.  8,  1.  6,  3.  2 and  6.  4 
seconds  and  formed  16  km  /sec  beams  in  twelve  directions,  evenly  spaced 
at  30  degree  intervals.  For  each  azimuth  and  integration  time  the  output  of 
each  detector  was  cpiantized  at  0.4  d!3  intervals  and  displayed  every  0.8 
seconds  in  the  same  form  as  used  for  the  long-period  data  in  Report  1. 

In  order  to  exclude  multiple  detections  from  the  same  event 
and  to  reduce  the  false  alarms  from  imperfectly  corrected  spikes,  a "dead 
time"  of  24  seconds  followed  the  time  at  which  th  ? detector  exceeded  each 
quantization  level.  During  this  time  no  more  detections  at  that  level  were 
allowed. 

The  variation  in  false  alarm  rate  from  one  azimuth  to  another 
was  small,  so  data  from  all  azimuths  at  a given  level  were  averaged  to- 
gether to  produce  the  false  .alarm  rates  shown  in  Figure  III  - 1 . Here  the 
number  of  detections  per  hour  is  plotted  versus  the  threshold  value  at 
which  those  detections  occurred.  Data  for  each  of  the  five  noise  samples 
studied  are  shown  for  the  Fisher  and  conventional  detectors  having  a 0.  8 
second  integration  time.  Since  the  24  second  dead  time  allowed  a maximum 
of  150  detections  in  one  hour,  the  curves  approach  this  value  at  tin-  lowest 
threshold  levels. 


Ill-  1 


nume 


I 


Two  broad  conclusions  may  be  drawn  from  these  curves. 

First,  the  conventional  detector  consistently  shows  a higher  alarm  rate 
than  does  the  Fisher  detector  at  the  same  threshold.  Second,  the  variation 
in  false  alarm  rate  from  day  to  day  is  greater  for  the  Fisher  detector  than 
for  the  conventional  detector.  Except  for  day  201,  the  conventional 
detector  is  substantially  the  same  from  one  day  to  another.  The  Fisher 
detector,  however,  can  vary  in  false  alarm  rate  by  at  leas,  a factor  of 
20  between  different  days. 

The  difference  in  alarm  rate  between  th^  detectors  is  rot 
significant  without  a knowledge  of  the  detection  probability,  and  cannot  be 
vised  to  conclude  that  one  detector  is  sviperior  to  the  other.  The  variation 
in  the  Fisher  detector's  false  alarm  rate  suggests  that  if  the  Fisher  detector 
was  superior  on  some  days  it  might  not  be  on  others,  depending  on  the  nature 
of  the  noise. 

The  root -mean  - sepia  red  amplitude,  averaged  over  all 
operating  channels  and  over  each  hour  of  noise  data  was  as  shown  in  Table 
III-  1 . Comparison  of  this  table  with  Figure  III- 1 shows  that  there  is  a general 
trend  for  days  with  high  rms  noise  levels  to  have  low  Fisher  detector  false 
alarm  rates;  the  results  on  day  201  are  the  only  exception  to  this  trend.  The 
large  difference  in  behavior  of  the  conventional  detector  on  day  20  1 from  its 
behavior  on  other  days  suggests  that  some  diife rent  noise  mechanism  was 
operating  on  that  day. 

In  Report  1 it  was  shown  that  spikes  in  the  data  reduce  the 
Fisher  detector  output.  When  numerous  small  spikes  or  large  amplitude 
excursions  have  been  present  on  noisy  days,  the  Fisher  detector  has  been 
observed  to  show  a lower  alarm  rate.  The  conventional  detector  averages 
these  spikes  ove  r all  channels  and  thus  is  relatively  insensitive  to  them.  Or 
day  20  1,  rapid  fluctuations  in  the  noise  power  about  its  mean  may  have  pro- 
duced the  observed  behavior  of  the  conventional  defector. 
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TABLE  III  - 1 

ROOT  MEAN  SQUARE  NOISE  AMPLITUDES 


Day 

RMS  Amplitude  (Digital  Counts) 

131 

4.  8 

141 

7.  5 

161 

4.  5 

183 

3.  8 

201 

8.  8 

HI -4 
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Figure  III-2  shows  false  alarm  rates  versus  level  for  the 
Fisher  detector  on  day  161,  for  a number  of  inleeration  times.  These 
behaviors  are  typical  in  that  the  spacinB  between  curves  is  maintained  from 
day  to  day,  although  the  curves  are  displaced  as  shewn  in  Figure  III-l. 

The  variability  of  the  noise  reported  here  means  that  a 
detector  using  a fixed  threshold  will  be  less  than  optimum.  For  example, 
if  a false  alarm  rate  of  5 per  hour  was  chosen  as  satisfactory,  experience 
on  day  131  would  indicate  that  a Fisher  detector  threshold  level  of  5,6 
dU  should  be  chosen.  Only  30  days  later,  on  day  161,  the  rate  at  this 
threshold  would  have  risen  to  about  17  per  hour,  an  unacceptable  figure. 

On  the  other  hand,  the  false  alarm  rate  on  day  141  would  have  been  about 
one  alarm  every  two  hours.  Such  a low  false  alarm  rate  is  of  course 
desirable,  but  on  this  day  a lower  threshold  would  have  realized  an  increased 

probability  of  event  detection  while  still  keeping  the  false  alarm  rate  within 
acceptable  bounds. 
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FISHER  AND  CONVENTIONAL  DETECTORS  FALSE  ALARM  RATES 
ON  DAY  161  AT  VARIOUS  INTEGRATION  TIMES 


SECTION  IV 

DETECTION  PROBABILITY 


In  this  section  expressions  for  comparing  and  evaluating 
automatic  detectors  are  presented , and  some  sample  results  are  shown. 

These  curves  relate  detection  probability  to  magnitude  for  a given  threshold. 
Since  the  alarm  rates  found  in  section  III  are  so  variable,  the  results  shown 
here  are  to  be  regarded  as  illustrative  of  the  method  and  to  form  a standard 
by  which  improved  versions  of  these  detectors  may  be  judged. 


All  the  events  described  in  Section  II-B  were  processed  with 
both  detectors.  The  largest  amplitude  of  the  detector  output  occurring 
within  36  seconds  of  the  predicted  arrival  time  on  three  beams  spaced  at 
15  degrees  and  centered  on  the  signal  azimuth  was  recorded.  The  magnitude 
assigned  to  the  event  was  that  listed  by  NORSAR  or  LASA.  Then  the  number 
of  detections  and  non-detections  at  each  magnitude  were  found  by  counting 
the  number  of  events  for  which  the  detector  output  did  or  did  not  exceed  an 
assumed  detector  level,  respectively. 

These  data  were  fitted  to  a model  adapted  from  one  developed 
by  Ringdal,  (1974).  His  model  assumes  that  the  probability  of  detecting 
an  event  of  magnitude  M is  given  by  the  cumulative  Gaussian  distribution 


M 


P(M)  = 


/■ 


- 

2 'o' 


dx 
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IV-1 


when  p and  (Tarn  parameters  to  be  fitted  to  tlu*  observed  detection  probability. 
Ibis  f 1 1 1 i n n is  done  so  as  to  maximize  tho  likelihood  fum  lion  associated  with 
equation  IV- 1.  The  function  P(M'  > s (ailed  the  detection  curve. 

This  theoi  , assumes  1 nat  the  probability  of  delecting  a 
sufficiently  small  event  is  zero.  For  the  automatic  detector  there  is  a 
probability  related  to  the  false  alarm  rate  that  the  detector  'evcl  will  exceed 
the  threshold  within  the  signal  arrival  time  window  due  to  noise.  This  pro- 
bability of  incorrectly  detec'ing  small  events  was  incorporated  into  the  maxi- 
mum likelihood  approac  h by  modifying  equation  IV- 1 to  read 


P(M)  K Ml 


IV  -2 


Here  K is  the  constant  detection  probability  at  small  magnitude,  and  is  found 
by  the  fitting  process  as  arep'  and  rr 1 . As  a practical  matter,  it  is  impossible 
to  estimate  ai  c irately  three  parameters  from  the  limited  data  sample 
available,  s<>  the  results  here  shouP  be  regarded  as  illustrations  of  the 
method. 

An  example*  of  a detection  curve  and  the  assocated  histogram 
of  detections  and  non-detections  are  shown  in  Figure  IV- 1 for  the  Fisher 
detector  at  0.  8 so'-oml  integration  time  and  threshold  level  of  5.  2 dB.  The 
small  magnitude  detection  orobability  is  .15,  and  MB50  and  SIGMA  arc  fi' 
and  (T  respectively.  Confidence  limits  for  this  curve  are  difficult  to  calculate 
when  variations  in  K art*  allowed,  and  those  shown  do  not  account  for  this 
variation.  This  problem  will  be  overt  time  in  the  future  by  fixing  K before 
maximizing  the  probability  function. 
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FIGURE  IV- la 


DISTRIBUTION  OF  PROCESSED  EVENTS  BY  m 


0.  8 SECOND  FISHER  DETECTOR  AT  LEVEL  5.  2 dB 


MAXIMUM  LIKELIHOOD  DETECTABILITY  CURVE 
0.8  SECOND  FISHER  DETECTOR  AT  LEVEL  5.2  dB 
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When  the  analyst  checks  the  claimed  detections  on  which  the 
curve  of  Figure  IV-lb  is  based,  he  will  reject  those  which  are  due  to 
noise  virtually  every  time.  Therefore,  this  curve  lies  above  the  true 
detection  curve  given  by 


P(M)  = / 2 a'  IV'3 

J 2 J e dx 

v2  no  - oo 

where  the  parameters  are  determined  from  the  model  equation  IV-2.  This 
curve  is  the  probability  of  confirmed  detections  as  a function  of  magnitude. 

It  and  equation  IV-2  are  shown  in  Figure  IV-2  for  the  data  of  Figure  IV-1. 

The  small  magnitude  intercept  of  the  observed  detection  curve 
is  the  probability  that  a detector  output  equal  to  or  greater  than  the  threshold 
occurred  at  least  once  in  three  24  second  intervals  on  three  adjacent  azimuths, 
due  to  noise  alone.  Let  the  probability  of  a false  alarm  be  P Per  24  seconds. 
Then  the  probability  of  no  false  alarm  over  all  three  azimuths  for  72  seconds 
is  (l-P^)11,  where  n = 3 if  the  data  are  perfectly  correlated  from  azimuth  to 
azimuth  and  = 9 if  they  are  perfectly  uncorrelated.  This  probability  of  no 
alarm  for  the  data  ot  Figure  1 V - 1 is 

1 - . 15  = . 85  = (1-P  .)n  r’-4 

L *T 

giving  a range  for  P of  0.  05  to  0.  02,  corresponding  to  a range  of  7.  5 to  3 false 
alarms  per  hour.  Reference  to  Figure  III  - 1 shows  that  this  rate  is  reasonable. 

The  uniformity  from  day  to  day  in  the  spacing  of  the  false 
alarm  curv  s for  various  integration  times  displayed  in  Figure  III  - 2 means 
that  a comparison  of  detectors  at  different  integration  times  can  be  made.  For 
this  purpose  the  detection  curve  for  the  Fisher  detector  using  a 6.4  second 
integration  time  was  calculated  assuming  a threshold  of  1.6  dB. 
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Detection  Probability 


FIGURE  JV-2 


MAXIMUM  LIKELIHOOD  DETECTABILITY  CURVE 
TRUE  AND  OBSERVED  FISHER  DETECTOR 
AT  LEVEL  5.  2 dB 
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This  level  should  give  a false  alarm  rate  comparable  to  that  achieved  with 
the  Fisher  detector  using  a threshold  of  5.  2 dB  and  an  integration  time  of 
0.  8 seconds,  according  to  the  discussion  of  Figure  III- 2-  The  detection 
histogram  and  curve  are  shown  in  Figure  IV- 3a  and  IV-3b. 

The  curve  in  Figure  IV-3b  is  indistinguishable  from  that  of 
Figure  IV-lb.  The  difference  in  small  magnitude  detection  probability,  0.  1 
in  Figure  IV-lb  and  0.  15  in  Figure  IV-3b,  is  no  larger  than  might  be 
expected,  given  the  small  number  of  events  in  the  data  sample.  Consequently, 
we  can  conclude  that  there  is  no  discernible  difference  between  the  Fisher 
detectors  using  6.4  and  0.8  second  integration  times  when  the  thresholds 
are  adjusted  to  give  a constant  false  alarm  rate.  A similar  experiment 
shows  that  there  is  likewise  no  difference  between  the  0.8  second  and  6.4 
second  conventional  detectors,  tor  levels  which  give  comparable  false  alarm 
rates. 
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Detection  Probability 


r * ■ 1 


DISTRIBUTION  OF  PROCESSED  EVENTS  BY  mb 
6.  4 SECOND  FISHER  DETECTOR  AT  LEVEL  1. 6 dB 


MAXIMUM  LIKELIHOOD  DETECTABILITY  CURVE 
6.  4 SECOND  FISHER  DETECTOR  AT  LEVEL  1.6  dB 
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CONCLUSIONS 


Fisher  and  conventional  beam  power  seismic  event  detectors 
similar  to  those  reported  on  previously  have  been  developed  for  short-period 
data  and  subjected  to  a preliminary  analysis  on  short-period  data  from  the 
Korean  Seismic  Research  Station. 

It  was  found  that  the  false  alarm  rates  produced  by  the  Fisher 
detector  when  using  fixed  detector  thresholds  varied  by  as  much  as  a factor 
of  20  from  day  to  day.  For  the  conventional  detector,  however,  variations 
in  the  false  alarm  rate  were  smaller  and  occurred  less  often.  The  Fisher 
detector  threshold  level  for  a given  false  alarm  rate  was  consistently  lower 
than  that  for  the  conventional  detector.  Increased  noise  power  generally  re- 
sulted in  lower  Fisher  detector  false  alarm  rates,  but  caused  little  change  in 
the  conventional  detector  alarm  rate.  Increasing  the  integration  time  of  the 
detectors  decreased  the  false  alarm  rate  in  a regular  manner. 

A total  of  185  events  w'hich  occurred  in  the  spring  of  1973  were 
analyzed  with  both  detectors.  Detection  curves  for  a few  fixed  threshold  levels 
ind  integration  times  were  developed  taking  into  account  the  possibility  of 
incorrect  detections  at  small  magnitudes.  These  curves  suggest  that  neither 
the  conventional  detector  nor  the  Fisher  detector  is  affected  by  changing  the 
integration  time  while  keeping  the  false  alarm  rate  constant. 

The  wide  variation  in  false  alarm  rates  at  a fixed  detector 
threshold  means  that  a detector  operating  with  a fixed  threshold  will  not  be 
optimum.  It  also  means  that  the  performance  of  such  Fisher  and  conventional 
detectors  cannot  be  compared  directly.  Optimum  detectors,  which  maintain 
fixed  alarm  rates,  could  be  used  for  this  purpose.  At  a constant  alarm  rate, 
whichever  detector  had  the  higher  probability  of  detection  would  unambiguously 
be  the  better  detector. 
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